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^Measurements  are  reported  of  the  growth  of  air  bubbles  in  water 
by  rectified  diffusion  at  22.1  kHz.  Values  of  the  threshold  acoustic 
pressure  amplitude  were  obtained  as  a  function  of  bubble  radii, s,  liquid 
surface  tension  and  gas  concentration.  Measurements  of  the  rati  of 
growth  of  bubbles  by  rectified  diffusion  were  also  obtained  as  a  func¬ 
tion  of  acoustic  pressure  amplitude  for  a  range  of  different  values 
of  the  liquid-vapor  surface  tension.  It  was  determined  that,  although  — = 
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both  the  threshold  and  the  growth  rate  were  in  agreement  with  theory  for 
nortuai  values  of  the  surface  tension  of  water,  the  addition  of  a  surfactant 
caused  the  observed  thresholds  and  growth  rates  to  deviate  from  the 
predicted  values.  Surface  wave  activity  that  could  increase  the  diffusion 
rate  by  acoustic  streaming  was  not  detected  at  low  radii  and  is  not 
thought  to  be  the  principal  mechanism  for  the  increased  diffusion.  Some 
possible  explanations  are  given  for  the  effect. 
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INTRODUCTION 


This  paper  concerns  Che  phenomenon  of  rectified  diffusion  of  gas 

into  air  bubbles  that  are  caused  to  pulsate  ip  a  liquid  by  the  action 

of  an  acoustic  field.  Free  bubbles  that  would  normally  dissolve  in 

water  may  be  made  to  grow  due  to  the  unequal  mass  transfer  across  the 

air-water  interface  during  bubble  oscillation.  The  phenomenon  seems 

to  have  originated  with  Harvey, ^  and  was  first  discussed  in  some  detail 
2 

by  Blake,  who  formulated  an  inadequate  theory  for  predicting  the  thresh- 

3  ^ 

old.  Attempts  by  Pode  and  Rosenberg  were  also  found  inadequate  to 

explain  some  of  the  early  measurements  by  Strasberg.”*  Hsleh  and  Ples- 

set,^  by  adding  a  convection  term  in  the  diffusion  equation,  considerably 

improved  upon  Blake's  approximate  results.  Their  theory  was  shown  by 
7 

Strasberg  to  be  in  good  agreement  with  his  measured  values.  Eller 

ft 

and  Flynn  ’  extended  the  theory  to  include  nonlinear  or  large  amplitude 
q 

effects  and  Safar  showed  that  the  Hsieh-Plesset  and  Eller-Flyum  treat¬ 
ments  were  essentially  equivalent  when  inertial  effects  were  not  neg¬ 
lected  in  the  Hsieh-Plesset  approach.  Eller  ’ ^  made  several  measure¬ 
ments  of  both  the  threshold  and  the  growth  rate  and  found  that  his  theory 
was  adequate  in  predicting  thresholds  but  was  unable  to  account  for  some 
very  large  growth  rates  that  he  observed.  He  suggested  that  the  rapid 

12 

rates  of  growth  may  be  accounted  for  by  acoustic  microstreaming.  Gould 
has  observed  bubbles  growing  by  rectified  diffusion  ^nd  has  found  agree¬ 
ment  with  Eller's  theory  in  the  absence  of  streaming.  By  observing  the 
bubble  directly  in  a  microscope  he  was  able  to  detect  the  onset  of  sur¬ 
face  oscillations  of  the  bubble  that  in  turn  introduced  significant 
acoustic  microstreaming.  This  streaming  greatly  enhanced  bubble  growth 


and  was  likely  the  explanation  lor  the  large  growth  rates  seen  by  Eller, 
13  14 

Attempts  by  Davidson  and  by  Kapustina  and  Statnikov  to  account  for 

12 

enhanced  growth  rates  by  miciostreaming  appear  to  be  inadequate. 

In  this  paper  measurements  will  be  presented  of  the  growth  of  air 
bubbles  by  rectified  diffusion,  for  various  values  of  the  surface  ten¬ 
sion  of  the  host  liquid,  the  acoustic  pressure  amplitude,  the  dissolved 
air  content  and  the  bubble  radius.  Comparisons  with  theoretical  pre¬ 
dictions  indicate  good  agreement  with  theory  for  both  the  threshold  and 
the  growth  rate  for  normal  values  of  the  surface  tension  of  water,  pro¬ 
vided  surface  wave  activity  of  the  bubble  is  avoided.  However,  the  ad¬ 
dition  of  a  surfactant  to  the  water  results  in  a  slight  reduction  in  the 
experimental  threshold,  compared  with  the  calculated  values,  and  in  a 
substantial  increase  in  the  rate  of  growth  of  the  bubble.  Although  the 
threshold  for  surface  wave  activity  is  lowered  for  reduced  values  of 
the  liquid  surface  tension,  an  explanation  for  the  increased  diffusion 
in  terms  of  bubble  surface  oscillations  appears  to  be  inconsistent  with 


the  observations. 


6 


1.  EXPERIMENTAL  APPARATUS  AA’D  PROCEDURE 


The  various  data  concerning  the  growth  ot  bubbles  by  rectified  dif¬ 
fusion  were  obtained  by  acoustically  levitating  the  air  bubbles  near  the 
antinode  of  an  acoustic  stationary  wave,  llus  technique  ’  , 

although  quite  common  for  many  applications,  has  been  considerably  re¬ 
fined  for  these  measurements.  A  diagram  of  the  experimental  apparatus 
is  shown  in  Figure  1.  The  stationary  wave  system  was  constructed  by 
cementing  a  hollow  glass  cylinder  between  two  matched  hollow  cylindri¬ 
cal  transducers,  fitted  with  a  flexible  pressure  release  diaphragm  on 
one  end  and  open  at  the  other.  The  composite  system  was  approximately 
?.5  cm  in  diameter  by  10  cm  in  height,  the  width  of  the  glass  in  the 
middle  about  2.5  cm.  This  system  was  driven  at  its  (r,G,z)  =  (2,0,2) 
mode  at  a  frequency  of  22.1  kHz. 

In  order  to  obtain  accurate  values  of  the  acoustic  pressure  ampli¬ 
tude  required  to  acoustically  levitate  a  bubble  at  a  specific  position, 
it  is  necessary  to  know  the  spatial  variation  of  the  stationary  acoustic 
wave  along  the  axis  of  the  cylinder.  This  variation  was  measured  with 
the  use  of  a  small  calibrated  probe  hydrophone  that  was  mounted  to  a 
micromanipulator.  It  has  been  determined  previously  that  stationary  wave 
systems  such  as  the  one  used  here  often  have  urstcrced  acoustic  pressure 
profiles  and  an  accurate  measurement  of  this  profile  is  necessary  for 
accurate  acoustic  pressure  amplitude  measurements.^  In  order  to  reduce 
effects  of  the  hydrophone  itself,  the  actual  acoustic  pressure  ampli¬ 
tude  measurements  were  made  by  a  small  pill  transducer  mounted  externally 
as  shown  in  Figure  1.  The  acoustic  pressure  amplitude  at  the  bubble's 
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position  was  determined  from  cne  voltage  output  of  the  pill  transducer 
and  the  known  spatial  variation  of  the  sound  field.  The  average  radius 
of  the  bubble  was  determined  by  first  measuring  its  terminal  rise  velo¬ 
city  through  the  host  liquid  and  then  by  iterating  the  following  equn 
17 

tion: 

+  ->.b  x  10~4  R^1'4}  (!) 

In  Eg.  1,  v  is  the  kinematic  viscosity  of  the  liquid,  u  is  the  terminal  j 

I 

velocity,  g  is  the  gravitational  acceleration  and  R^  =  2  I^u/v  is  the  j 

Reynolds  number.  Two  refinements  were  made  to  inctcase  the  accuracy  j 

of  the  radius  measurements.  First,  several  fiducial  lines  were  intro¬ 
duced  into  the  cathetometer  so  that  rise-times  could  be  taken  over  a 
variety  of  distances  and  thus  would  not  be  excessively  short  or  long. 

Second,  the  rise-time  measurements  for  an  individual  bubble  was  intro¬ 
duced  directly  into  a  computer  program  that  corrected  foi  the  dissolution 
of  the  bubble  when  the  sound  field  was  off  and  the  terminal  velocity 
measurement  was  being  made.  For  small  bubbles,  say  20  pm  or  less,  an 
interval  of  a  few  seconds  used  to  measure  the  bubble's  rise  velocity 
could  result  In  errors  of  10%  or  more  in  the  bubble  radius  due  to  bub¬ 
ble  dissolution. 

The  surface,  tension  of  the  host  liquid  was  varied  by  addition  to 

distilled  water  of  small  amounts  of  a  commercially  available  surfactant, 

ethexylated  octaphenol,  marketed  under  the  brand  name  of  Photoflo,  and 

commonly  used  in  photography.  The  surface  tension  was  measured  with  a 

18 

Du  Nouy  ring  tensiometer 

It  was  discovered  that  small  variations  in  the  temperature  of  the 
host  liquid  could  result  in  substantial  changes  in  the  rectified  diffusion 


R  = 


9vu 

2g 
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1  +  0.20  R 


0.63 


thicshold  and  growth  rate  ol  the  bubble.  As  will  be  shown  in  Section 
111,  a  small  change  in  temperature  could  result  in  a  signil leant  ei- 
fect  due  to  the  strong  variation  in  the  dissolved  gas  concentration  with 
temperature.  Accordingly,  it  was  necessary  to  surround  the  test  cell 
with  a  constant  temperature  bath,  see  Figure  1,  to  insure  temperature 
stability,  and  equilibrium  gas  concentration.  The  temperature  of  the 
host  liquid  was  constantly  monitored  to  insure  no  change  in  temperature 
during  the  measurements,  or  just  as  importantly,  no  change  from  the  tem¬ 
perature  at  which  the  liquid  was  brought  to  equilibrium  gas  concentra- 


In  general  it  is  estimated  that  the  absolute  accuracy  of  the  acous¬ 
tic  pressure  amplitude  measurements  was  good  to  within  5%,  and  consistent 
to  within  2%.  The  radius  measurements  were  slightly  better  except  tor 
small  radii,  <  ?Sum.  Surface  tension  measurements  tended  to  be  subjected 
to  contamination  during  an  experiment,  but  were  consistent  to  within  3% 
and  of.  similar  absolute  accuracy.  Dissolved  gas  concentration  was  not 
measured  directly,  but  considerable  care  was  made  to  bring  the  liquid  to 
equilibrium  by  vigorous  agitation  of  the  host  liquid,  for  several  minutes, 
and  then  maintaining  the  temperature  to  within  0.2°C. 

Bubbles  were  introduced  into  the  test  cell  by  increasing  the  gain  on 


the  amplifier  until  the  t.  reshold  for  gaseous  cavitation  was  achieved. 

By  suitable  manipulation  of  the  acoustic  pressure  amplitude,  bubbles  of  the 
desired  size  could  be  acquired  within  a  few  minutes. 


-  v.  * 


HYDROPHONE 


FJG  i.  Block  diagram  and  sketcd  of  experimental  arrangement 


11.  EQUATIONS  FOR  BUBBLE  GROWTH 


In  this  section  we  obtain  the  equations  that  describe  the  change  in 

size  of  an  air  bubble  that  is  present  in  a  liquid  and  experiencing  a 

periodic  fluctuation  in  pressure  described  by  P.  cos  ul .  where  P  is  the 

A  t\ 

aeoustic  pressure  amplitude  and  to  is  the  angular  frequency  of  the  oscillation. 
The  bubble  will  either  grow  or  dissolve  depending  upon  a  variety  of 
conditions;  we  examine  the  conditions  necessary  for  size  stabilization 
(threshold  for  rectified  diffusion)  and  for  change  in  size  of  the  bubble 
as  a  function  of  time  (growth  rate). 

Previously,  expressions  have  been  obtained  for  the  threshold  and  the 
growth  rate  by  Hsich  and  Plesset6,  Strasberg7,  Saiar°  and  Eller10  for 
isothermal  pulsations  ot  the  bubble,  and  by  Eller11  for  adiabatic 
pulsations.  Our  experimental  measurements  at  22.1  kilz  are  for  bubbles 
that  span  the  gap  between  purely  isothermal  pulsations  (R^  <  30pm)  and 
those  that  approximate  adiabatic  pulsations  (RQ  >  70pm).  Accordingly,  we 
have  attempted  to  obtain  expressions  for  the  growth  rate  and  the  threshold 
acoustic  pressure  amplitude  that  apply  to  the  entire  range  of  bubble 
sizes . 

11 

Eller  has  simplified  the  more  general  theory  and  has  obtained  an 
equation  for  the  rate  of  change  in  the  number  of  moles  of  gas  in  a  bubble 
to  be 

-  4n  DR  <  R/R  >  C  H  (2) 

dt  o  o  o 

where  a  time-dependent  term  has  been  neglected,  P  is  the  diffusion  constant, 

R,  are  the  instantaneous  and  equilibrium  values  of  the  bubble  radius, 

C0  is  the  "equilibrium"  or  "saturation"  concentration  of  the  gas  in  moles/unit 


volume , 
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where  C.  is  the  concentration  of  dissolved 

a 
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where  o  is  the  surface  tension  of  the  host  fluid  and  n  is  the  polytropic 
exponent.  It  is  further  assumed  that  the  bubble  radius  obeys  the  equation 
of  motion 

RA  +  f  R2  +  p"1  [I2-  -  (P  +  I2)  (i£)3'1+  P  +  P.  cos  wt]  ■  0  (5) 

*.  K  O  K  K  O  A 

0 

where  p  is  tlie  density  of  the  liquid  and  the  dot  above  the  symbol  implies 
differentiation  "ith  respect  to  time.  Following  Eller"*"*  we  seek  an 
approximate  solution  of  Eq.  5  by  an  expansion  in  (P./P  ),  namely, 

A  O 
P 

R/R  *  1  -  a(~)  cos  mt  +  a2  K  (P./P  )2  +  ...  (6) 

o  P  A  o 

o 


whore  we  have  determined 


a 


K. 


[3n  (1  +  4o/3r  R  -  62)]'1 
o  o 

(3n  +  1  -  p2)/4  +  4o/3P  R 
_  _ o  o 

1  T  3P  R 
o  o 


(7) 

(3) 
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and 


2  2 

..  pui  R 

6  »  -TT*- 
3nP 

o 


(9) 


The  time  averages  ^R/Ro^  ,  ^(R/P.o)A^and  ^(R/Rq) A  (P^/P^)  ^will  be  required 
later  and  are  given  by 


A/R  >  =  1  +  K  q2(P/P  )2  ,  (10) 

O  A  O 

<"(R/R  )A>  =  1  +  (3  +  4K)  u2  (PA/P  )2  ,  (11a) 

'  Of  A  o 

and  /(R/R  )A  (P  /P  )\  =  [l+(4-3n)Ka2(P,/P  32](l+2o/R  P  )  .  (lib) 

*  c  go*  A  o’  oo 

We  can  relate  the  equilibrium  value  of  the  radius  to  the  number  of  moles  of 
gas  by  the  equation 

F  +  ~  =  3nP.  T/4tiR  3  ,  (12) 

o  R  g  o 

o 

where  R^  is  the  universal  gas  constant  and  T  is  the  absolute  temperature. 
Combining  Eq.  12  with  Eq.  2  and  Eq.  3,  we  obtain  an  expression  for  the 
rate  of  change  of  the  bubble  radius  with  time 


dR 
_ o 

dt 


~  {<R/rN  (1+W3P  R  )-i 

K  »  O »  O  O 

o 


<(R/R  >  (P  /P  )  > 


(13) 


where  d  =  R  TC  /P  . 

goo 

In  Eq.  13,  <R/R3>,  <^(R/Ro)A^  and  ^(R/R^)A (P  /P^)^  are  given  by 
Eqs .  10  and  11  respectively. 

The  threshold  for  rectified  diffusic  itained  by  setting  dR^/dt  -■  0 

in  Eq .  13,  and  is  given  by 


■J*W,S*«*,*  wbtttitn'..  .. 
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P  1  +  2a/R  1’  -  C  JC 

o  t _  o  o  i  o 

^  l  V 


1/2 


a  •(3+AK'IC  /C  -  C1+20/R  P  )(4-3n)KJ 
x  o  oo 


13 


(14) 


The  polytropic  exponent  n  is  given  by  a  set  of  complicated  expressions 

19  20  21 

by  Eller  based  on  the  work  of  Devin  ,  and  by  Prosperetti  ,  using  a  some¬ 
what  different  technique.  For  completeness  we  list  the  expression  obtained 
by  Eller^ 


/,,j  »-l  r ,  .  3(7-1)  ,sinh  X  -  sin  XN1“1 
n  =  7(l+dt)  [1  +  —  (coihlc  ~o7x)1 


(15) 


where 


,  , ..  fX(sinh  X  +  sin  X)  -  2(cosh  X  -  cos  X)  , 

d  «  3(7-1)  - - - - - - - - ],  (16) 

X  (cosh  X  -  cos  X)  +  3(y-l)X  (sinh  X  -  sin  X) 


r 


T 

I 

i. 


and  X  =  R  (2u)/DJ 
o  1 


(17) 


In  these  equations  7  is  the  ratio  of  specific  heats  and  D,  =  K , /P,C  ,  where 

1  1  1  Pj 

is  the  thermal  conductivity  of  the  gas  ir.  the  bubble,  is  he  density 
of  gas  and  Cp  ,  is  the  specific  heat  at  constant  pressure  for  the  gas. 

If  the  polytxopic  exponent  n  is  set  equal  to  1.0  in  Eqs.  7-9,  then 
the  expressions  for  the  growth  rate,  Eq.  13,  and  for  the  threshold,  Eq.  14, 
reduce  to  those  derived  for  the  isothermal  case  by  Eller. ^  If  n  is  set 
equal  to  7,  then  Eqs.  13  and  14  reduce  approximately  to  those  obtained  by 
Eller^  for  the  adiabatic  case.  Eller^  neglected  the  eff«e.t  of  surface  tension 
and  s<’t  37  =  4.0.  approximations  not  made  in  this  study.  Comparisons  of  the 
equations  obtained  for  the  threshold  by  various  authors,  together  with  some 
experimental  measurements,  will  be  given  in  section  III. 


The  constants  that  were  used  in  the  above  equations  and  applicable  to 

this  experiment  are  f  =  22.1  kHz,  Y  =  1.4,  =  0.20  cm^/sec,  Pq  =  1.01  X 

2  3 

10  dyn/cm  ,  and  3  =  1.0  g/cm  .  The  constants  D  and  d  were  corrected  for 


temperature  variation;;  by  the  following  equations:  D  =  (6.15  X  10  T  - 

15.6  X  10  ">)  cm^/sec  and  d  =  1.02  X  10  ^  -  2.80  X  10  ^T,  T  being  the 

absolute  temperature  in  degrees  Kelvin.  The  data  used  in  the  equations 

for  the  variation  of  D  and  d  with  temperature  were  taken  from  the  Inter 
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national  Critical  Tables 


III.  EXPERIMENTAL  RESULTS 


Measurements  of  the  threshold  and  the  rate  of  growth  of  bubbles  by 
rectified  diffusion  were  made  for  a  variety  of  experimental  conditions.  In 
this  section  we  examine  the  results  of  the  measurements  and  the  comparisons 
with  theory. 

A.  Experiments  at  normal  surface  tension 
1.  Threshold  measurements 

Figure  2  shows  measurements  of  the  threshold  for  rectified  diffusion 

of  air  bubbles  in  water  at  a  surface  tension  of  68  dyn/cm  (distilled  water 

having  an  initial  surface  tension  of  72  dyn/cm  allowed  to  stand  in  our 

laboratory  at  a  temperature  of  25°C  for  a  few  hours  would  reach  an 

equilibrium  value  of  68  dyn/cm  no  matter  what  precautions  we  would  take 

against  contamination).  The  values  of  the  threshold  acoustic  pressure 

were  obtained  by  levitating  an  air  bubble  and  monitoring  its  radius  as 

a  function  of  time.  An  acoustic  pressure  amplitude  that  resulted  in  a 

radius  change  of  less  than  1%  over  a  time  interval  of  a  few  minutes  was 

considered  to  be  a  threshold  value.  Since  this  value  is  a  position  of 

metastable  equilibrium,  considerable  patience  is  required  in  "tracking 

down"  a  point.  Also  shown  in  Fig.  2  are  a  set  of  curves  represent ing  the 

various  theories  that  can  be.  applied  to  the  measurements.  Eller^  and 
9 

Safar  have  developed  equations  that  assume  isothermal  pulsations  of  the 
bubble;  Eller^  has  also  extended  his  work  to  include  adiabatic  pulsations. 

It  is  seen  that  Eq.  14 

more  closely  approximates  the  data  than  the  two  extreme  treatments,  and 


BUBBLE  RADIUS  (Mm) 

FIG.  2  Values  of  the  threshold  acoustic  pressure  amplitude  for  rectified  diffusion  as  a  function 
of  radius  for  a  surface  tension  of  68  dyn/cra  and  a  frequence  of  22.1  kHz.  The  dashed 
curves  indicate  the-  calculated  thresholds  for  the  limiting  cases  of  isothermal  and  adiaheti- 
pulsations.  Th.:  -.'id  curve  is  calculated  from  Eq .  14. 


17 


that  good  agreement  with  theory  is  obtained  for  the  threshold  for  rectified 
diffusion  for  air  bubbles  in  pure  water. 

Due  to  an  initial  inconsistency  in  our  measurements  we  examined  the 
effect  of  slight  changes  in  temperature  on  the  threshold.  For  the 
measurements  shown  in  Fig.  3,  the  following  procedure  was  employed.  A 
sample  of  water  was  vigorously  agitated,  brought  to  saturation  at  a 
particular  temperature,  and  then  placed  into  the  test  cell  within  the 
constant  temperature  bath,  set  at  the  same  temperature.  After  a  settling 
period  of  an  hour  or  so,  the  bath  was  raised  or  lowered  several  degrees 
in  temperature,  until  the  host  liquid  temperature  changed  a  few  degrees, 
and  then  the  bath  temperature  was  manipulated  to  maintain  this  new  temp¬ 
erature.  Measurements  of  the  threshold  were  then  made  at  this  temperature. 
Using  handbook  values  of  the  dissolved  gas  concentration  as  a  function  of 
temperature,  and  assuming  no  gas  diffusion  into  or  out  of  the  liquid  during 
this  period,  values  of  the  gas  oversaturation  or  undersaturation  of  the 
liquid  could  be  obtained.  It  is  seen  in  Fig.  3  that  a  rapid  temperature 
change  of  1°C  could  be  easily  detected,  and  temperature  changes  of  2.3°C, 
resulting  in  a  gas  oversaturation  of  approximately  4%,  could  result,  in  s 
dramatic  change  in  the  threshold.  Comparisons  with  theory  are  good 
considering  the  inexactness  of  the  technique.  It  is  seen  that  if  special 
care  is  not  taken  to  maintain  the  gas  content  of  the  liquid  at  saturation, 
significant  errors  can  be  introduced  into  the  measurements  by  a  slight 
drift,  in  temperature. 

2.  Growth  rate  measurements 

Measurements  were  also  made  of  the  rate  of  growth  of  bubbles  by 
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rectified  diffusion.  Figures  4  and  5  show  measurements  of  the  radii  of 

bubbles  as  a  function  of  time.  In  Fig.  4,  the  radius-time  curves  of  four 

bubbles  are  shown,  each  with  an  initial  radius  of  35  microns,  but  for 

different  values  of  the  acoustic  pressure  amplitude.  In  Fig.  5,  the 

radius-time  curves  of  three  bubbles  are  shown,  each  with  the  same  value 

of  the  acoustic  pressure  amplitude,  but  with  different  values  of  the  initial 

radius  of  the  bubble.  The  concept  of  a  threshold  pressure,  cr  a  threshold 

radius  is  readily  seen  from  the  figures.  Pitted  on  the  graphs  are  theoretical 

curves  generated  by  numerically  integrating  Eq.  13.  It  is  seen  that  good 

agreement  is  obtained  with  theory  for  the  growth  rates  of  the  bubbles 

provided  surface  waves  are  not  obviously  present.  In  Fig.  4,  for  the 

bubble  with  the  largest  acoustic  pressure  amplitude,  a  rapid  change  in 

the  growth  rate  was  observed  near  a  radius  of  60  microns.  It  was  also 

observed  that  the  scattered  light  from  this  bubble  changed  character, 

becoming  more  intense  and  beg.  n  fluctuating.  The  bubble  soon  began  large 

23 

scale  oscillations  that  have  previously  been  called  "dancing  motion." 

V/e  have  designated  this  point  on  Fig.  4  by  an  arrow  and  the  letters  SW 
to  Indicate  our  belief  that  the  bubble  was  undergoing  surface  wave 
oscillations. 

Figure  6  shows  the  results  of  a  study  of  the  growth  rate  as  a  function 
of  acoustic  pressure  amplitude  for  a  given  value  of  the  bubble  radius. 

The  experimental  points  in  Fig.  6  were  obtained  by  measuring  the  slope 
of  the  radius-time  curves  of  bubbles  at  a  radius  of  45  microns.  For 
example,  the  top  set  of  data  in  Fig.  show  a  bubble  that  has  a  growth 
rate  of  4.3  microns  per  100  seconds  at  a  radius  of  45  microns.  The 
data  in  Fig.  6  also  show  negative  growth  rates,  that  is,  for  bubbles 
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TIME  (SEC) 

Variation  of  the  bubble  radius  with  time  for  a  liquid  surface  tension  of  68 
dyn/cra,  and  for  various  values  of  the  acoustic  pressure  amplitude.  The  arrow 
indicates  the  detection  of  obvious  surface  waves.  The  liquid  was  saturated  v 
gas  and  the  solid  lines  are  calculated  by  numerical  integration  of  Eq.  13 


Variation  of  the  bubble  radius  with  time  for  a  liquid  surface  tension  of 
68  dyn/cm  and  an  acoustic  pressure  amplitude  of  0.27  bar.  The  liquid  was 
saturated  with  gas.  The  solid  lines  are  calculated  by  numerical  integration 


Variation  in  the  growth  rate  of  air  bubbles  by  rectified  diffusion  with 
acoustic  pressure  amplitude  for  a  surface  tension  of  68  dyn/ctn,  a  gas 
saturated  liquid,  and  for  a  radius  of  kb  pm.  The  solid  line  is  calculated 
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that  were  below  threshold  and  were  dissolving.  The  smooth  curve  in  Fig. 

6  i6  calculated  iron  Eq.  13  with  Rq  “  45  microns,  and  shows  good  agreement 
with  the  experimental  values. 

It  is  thus  seen  from  the  above  measurements  that  so  long  as  there 
is  no  surface  waves  present  on  the  bubble,  both  the  threshold  and  the 
growth  rate  of  bubbles  by  rectified  diffusion  can  be  adequately  predicted 
by  theory  for  the  modest  range  of  variables  studied  here. 

B.  Experiments  at  reduced  surface  tension 

Measurements  were  made  of  the  threshold  and  the  growth  rate  of  bubbles 
in  water  for  which  the  surface  tension  was  reduced  by  the  addition  of 
small  amounts  of  a  surface  active  agent  or  surfactant.  These  measurements 
tended  to  give  lack  of  agreement  with  theory  for  both  thresholds  and 
^H'wui  rates . 

1 .  Threshold  measurements 

Measurements  similiar  to  those  shown  in  Fig.  2  were  made  for  a 

variety  of  liquid  surface  tensions.  The  procedure  for  the  measurements 

was  identical  to  those  for  pure  water,  except  that  small  amounts  of 

surfactant  was  added  to  the  liquid  after  it  had  been  placed  into  the 

test  cell  and  had  reached  temperature  equilibrium.  A  few  drops  of 

3 

surfactant  added  to  the  200  cm  volume  of  water  in  the  levitation  trans¬ 
ducer  was  sufficient  to  reduce  the  surface  tension  from  68  dyn/ern  to 
32  dyn/cm. 

Figure  7  shows  the  results  of  measurements  of  the  rectified  diffusion 
threshold  as  a  function  of  liquid  surface  tension  for  two  values  of  the 
bubble  radius.  Also  shown  are  the  theoretical  predictions  from  Eq.  14. 
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It  is  seen  that  the  data  are  approximately  represented  by  the  theoretical 
curves  but  lower  values  of  the  surface  tension  tend  to  give  values  of 
the  threshold  acoustic  pressure  amplitude  that  are  consistently  less  than 
theory.  As  an  example,  the  lower  6  data  points  on  the  lower  curve  in 
Fig.  7  average  to  a  threshold  of  approximately  0.12  bar,  with  a  standard 
deviation  of  0.02.  The  theoretical  value  for  the  threshold  at  a  surface 
tension  of  37  dyn/cm  is  0.165  bar.  Thus,  the  six  experimental  points 
average  to  a  value  that  is  slightly  more  than  3  standard  deviations  from 
the  predicted  value.  At  best  it  can  he  said  that  th’jre  are  indications 
that  the  threshold  is  lowered  with  addition  of  a  surfactant  to  the  host 
liquid. 

2 .  Growth  rate  measurements 

Measurements  were  also  made  of  the  growth  rate  of  bubbles  by 
rectified  diffusion  for  various  values  of  the  liquid  surface  tension. 
Figure  8  shows  measurements  of  the  radius-time  curves  of  three  bubbles, 
each  with  an  initial  radius  of  35  microns,  but  having  different  values 
of  the  acoustic  pressure  amplitude.  The  measured  surface  tension  for 
these  data  was  32  dyn/cm.  Plotted  on  the  graph  are  the  theoretical 
curves,  taken  from  numerical  integration  of  Eq.  13,  that  predict  the 
radii  of  the  bubbles  as  a  function  of  time.  Note  that  the  measured 
growth  rates  are  significantly  larger,  for  those  bubbles  that  are  growing, 
than  the  predicted  values. 

In  Figure  9  we  have  collected  the  growth  rate  measurements  for  three 
different  values  of  the  liquid  surface,  tension  as  a  function  of  acoustic 
pressure  amplitude.  Plotted  on  the  graph  are  the  calculated  values  from 
Eq.  13.  It  is  seen  that  substantial  increases  in  the  growth  rate  from 


that  predicted  by  theory  are  observed  when  a  surfactant  in  added  to  the 
water. 

In  Figure  10,  the  measured  growth  rate  for  a  given  acoustic  pressure 
amplitude  and  for  a  given  radiue  is  shown  as  a  function  of  surface 
tension.  Substantial  disagreement  is  observed  between  the  experimental 
points  and  the  predicted  values,  shown  as  the  curve,  taken  from  Eq.  13. 


SURFACE  TENSION  (DYN/CM) 

FIG.  10  Variation  in  the  growth  rate  of  air  bubbles  by  rectified  diffusion  with  surface 

tension  for  a  radius  of  50  um,  a  gas  saturated  liquid,  and  an  acoustic  pressure 
amplitude  of  0.22  bar.  The  solid  line  is  calculated  from  Eq.  13. 


IV.  DISCUSSION 


Eller^’^  had  found  that  observed  growth  rates  were  larger  than 

predicted  even  though  observed  and  predicted  thresholds  were  about  the 

same.  He  suggested  that  a  possible  explanation  for  the  increased  growth 

rate  was  due  to  acoustic  microstreaming  that  greatly  Increased  the  avail- 
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ability  of  gas  for  diffusion  into  the  bubble.  Gould  examined  the  e. ' uct 
of  microstreaming  on  diffusion  rates  and  concluded  that  the  onset  of 
surface  wave  activity  did  indeed  induce  streaming  which  markedly  increased 
the  diffusion  rate.  Further,  he  found  essential  agreement  for  the  pre¬ 
dicted  and  observed  growth  rates  when  microstreaming  was  absent.  The 
results  of  this  experiment  corroborate  these  findings.  However,  it  has 
been  discovered  in  this  experiment  that  at  reduced  surface  tension, 
observed  and  predicted  growth  rates  began  to  increasingly  diveige  as  the 
surface  tension  is  reduced.  It  was  first  thought  that  surface  wave 
activity  was  responsible  for  the  increased  growth  rates  observed.  However, 

]  2 

in  the  observations  by  Gould  (see  for  example  his  Figs.  7  and  8)  these 

surface  oscillations  were  accompanied  by  such  a  rapid  increase  in  the 

growth  rate  that  a  readily  observable  kink  occurred  in  the  radius-time 

curve.  Figures  4  and  8  of  this  study,  which  are  similar  to  over  100 

other  radius-time  observations  made,  show  cases  of  similar  kinks  in  the 

growth  curves,  but  only  after  the  observed  growth  rate  has  deviated 

significantly  from  theory  for  the  lower  surface  tensions. 

It  is  possible,  to  make  an  estimate  of  the  threshold  acoustic  pressure 

23 

amplitude  required  for  generating  surface  waves.  Eller  and  Crum  , 

6  24 

Hsieh  and  Plesset j  and  Benjamin  have  made  calculations  of  the  thres¬ 
hold  for  surface  wave  production.  An  approximate  expression  was  given 
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by  Hsieh  and  Plesset  that  agrees  reasonably  well,  at  a  radius  of  45 

23 

microns,  with  the  more  detailed  calculations  of  Eller  and  Crura  .  The 
threshold  acoustic  pressure  amplitude  required  to  generate  n  =  2 
oscillations  is  given  by 
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For  a  radius  of  45un  and  a  surface  tension  of-t>8  dyn/crc,  P  =  0.49  bar: 

A 

for  a  surface  tension  of  58  dyn/cm,  P^  =  0.42  bar;  and  for  a  surface 

tension  of  32  dyn/cm,  P^  ~  0.23  bar.  It  is  seen  from  Fig,  9  that  the 

deviation  of  the  measured  growth  rate  from  the  predicted  value  occurs 

before  these  surface  wave  thresholds  are  reached.  It  is  noted,  of  course, 

that  the  surface  wave  thresholds  listed  here  are  calculated  values,  and 

actual  thresholds  may  be  somewhat  less.  Some  measured  values  of  the 
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threshold  for  surface  wave  generation  do  exist  '  ,  but  there  is  con¬ 

siderable  scatter  in  the  results  by  Gould^2,  and  no  mode  identification 

23 

in  the  results  of  Eller  and  Crum.  Precise  measurements  of  the  thres¬ 
hold  for  surface  wave  generation  are  required  in  order  to  reject  or  accept 
this  explanation  tor  the  increased  diffusion  at  low  surface  tensions. 

A  second  possible  but  rather  speculative  explanation  is  now  proposed 
to  explain  the  increased  diffusion. 

It  has  been  known  for  some  time2’’  ’ 22  that  surface  active  agents 
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can  reduce  the  evaporation  rate  of  water.  Mansfield  has  shown  that 

the  water  loss  by  evaporation  from  large  lakes  may  be  reduced  by  as  much 

as  50%  by  addition  of  an  appropriate  monomolecular  film.  Archer  and 
28 

La  Mer  found  that  the  rate  of  evaporation  may  be  reduced  by  a  factor 
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of  10  by  the  application  of  a  fatty  acid  monolayer.  Jn  an  attempt  to 

29 

account  for  these  large  reductions  in  the  evaporation  rate,  Barnes  and 

co-workers  calculated  the  evaporation  resistance  assuming  that  water 

molecules  could  only  evaporate  through  holes  in  the  monolayer  film. 

With  a  reasonable  estimate  of  the  ratio  of  hole  area  to  monolayer  area 

they  were  able  to  calculate  values  of  the  evaporation  resistance  that 

26 

approximated  the.  measured  values.  Mansfield  also  found  that  expansion 
of  the  file,  by  as  little  as  10%  reduced  the  evaporation  resistance 
considerably,  while  compression  of  the  film  had  little  effect  on  the 
evaporation  resistance. 

The  existence  of  surface  active  materials  present  on  air  bubble 

30  31 

surfaces  have  been  previously  suggested  as  affecting  the  dissolution,  ’ 
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and  the  resonance  frequency  of  Small  air  bubbles.  Elder  visually 

observed  the  presence  of  a  layer  of  surface  active  materials  on  air 

34 

bubbles  and  its  effect  on  acoustic  microstreaming.  Liebermann  examined 

the  effect  of  surface  contamination  in  reducing  the  diffusion  of  a 

stationary  air  bubble  that  was  allowed  to  slowly  dissolve  in  water. 

However,  he  states  that  he  was  unable  to  find  a  significant  (by  a  factor 

of  two)  eEfect  on  the  diffusion  rate  by  adding  surface-active  agents  such 

as  detergents.  He  did  find  that  air  bubbles  left  a  solid  particulate 

residue  when  dissolved  on  a  glass  slide. 

It  appears  that  some  facts  point  toward  the  presence  of  a  surface 

active  monolayer  as  a  possible  explanation  for  the  large  growth  rates 

26 

observed  in  rectified  diffusion.  The  observation  by  Mansfield  that 


expansion  of  a  film  allows  more  evaporation  than  when  the  film  is  compressed 
is  directly  applicable  to  the  case  of  a  pulsating  air  bubble.  If  the  film 


allows  more  pend  ration  of  molecules  on  expansion  than  on  compression,  a 

rectification  oi  diffused  air  would  occur,  in  agreement  with  the  above 

observations  on  rectified  diffusion.  Further,  if,  on  expansion,  the 

film  is  broken  up,  the  time  required  for  the  film  to  repair  itself  is 
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long  compared  to  the  period  of  the  acoustic  cycle.  Lord  Rayleigh, 

in  an  ingenious  experiment  performed  in  1890,  showed  that  surface  active 

agents  needed  at  least  10  msec  for  their  migration  to  the  surface  and 

their  resultant  reduction  in  surface  tension.  On  the  other  hand,  the 

compression  may  cause  the  repair  mechanically  so  that  no  diffusion  of  the 

film  from  the  liquid  substrate  is  required.  At  any  rate,  it  seems  logical 

to  assume  that  if  evaporation  rates  can  be  retarded  by  four  orders  of 

magnitude,  diffusion  rates  can  be  affected  enough  each  cycle  to  cause 

a  significant  increase  in  the  amount  of  gas  contained  within  the  bubble. 

34 

The  observations  by  Liebermann  that  diffusion  rates  may  be  slightly 

affected  by  surface  contaminants  support  this  view.  With  the  rapid 

pulsation,  22.1  kHz,  very  little  rectification  is  required  each  cycle 

to  significantly  affect  the  growth  rate.  Further,  if  the  film  has  a 

rectifying  effect,  it  should  lower  the  threshold  for  rectified  diffusion. 

Although  the  measurements  do  not  show  an  obvious  and  distinct  lowering, 

there  are  indications  of  a  threshold  reduction  in  these  measurements. 

Finally,  the  absence  of  any  kinks  in  the  radius-time  curves  imply  no 

immediate  change  in  boundary  conditions,  as  would  occur  if  a  bubble 

underwent  surface  wave  oscillations. 

A  third  explanation,  and  perhaps  the  most  plausible  one,  is  that 

some  raicrostreaming  can  be  induced  without  surface  wave  generation.  It 

33 

was  observed  by  Elder  ,  for  example,  in  his  studies  of  microstreaming 
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that  when  surface  active  agents  were  added  to  water,  a  thin  film  would 
form  which  tended  to  present  a  no-slip  boundary  condition  tor  the  bubble. 
This  boundary  layer  effect  increased  the  microstreaming  until  it  was 
broken  up  by  the  bubble  pulsations.  It  is  possible  that  the  addition  of 
small  amounts  of  surfactant  induces  microstreaming  without  surface 
waves.  This  streaming  would  increase  the  growth  rate  for  all  values  of 
the  acoustic  pressure  amplitude  and  bubble  radius,  and  would  show  no 
distinct  inception  threshold. 
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